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TRAPPED-PARTICLE RADIATION ENCOUNTERED BY 


SELECTED LOW-ALTITUDE SKYLAB MISSIONS: 

SPACE STATIONS IN AN AMBIENT 
RADIATION ENVIRONMENT 


Foreword: 

At the request of The Office of Manned Space Flight/Code MT, and in order to 
supply the Office with necessary data and information that would enable decisions 
to be made about the "Space Platform" program during the early planning stages 
of the mission, a special radiation study was conducted for a series of selected 
trajectories. 

This report contains the results of our study; it presents an analysis of the data 
obtained with interpretive comments and evaluations, including all pertinent 
background information and explanations. 


Trajectory Generation and Conversion: 

For the specified flight-paths, orbit tapes were generated from the latest 
ephemerides, taking into account normal perturbations (drag, radiation pressure, 
planetary and lunar perturbations) wherever applicable, with an integration step- 
size of one minute, in all cases. 

Each trajectory was generated for a sufficiently long flight-time so as to insure 
an adequate sampling pf the ambient environment; on account of their period, 
which determines the orbit-orecession, the following circular flight paths of 
48-hour duration were produced: 
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The orbits were subsequently converted from geocentric polar into magnetic 
B/L coordinates with Mcllwain's INVAR program of 1965, utilizing the GSFC 
(12,/1966) geomagnetic field model by Cain et al., calculated for the epoch 
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1971.0, where B is the field strength and L is the geocentric distance to the 
intersect of the field line with the geomagnetic equator. 

A circular trajectory with high inclination (i > 55°) traverses the terrestrial 
radiation belts twice during each revolution. The vehicle thus executes a 
transverse motion in L-space, passing successively through a region of low 
L- values (1.0 < L £ 2.0) and of high L-values (2.0 £ L £ 6.6), commonly 
referred to as the inner zone and the outer zone. 


Models of the Trapped Particle Environment and Temporal Variations : 

Orbital flux integrations were performed with Vette's current environment models, 
the AE2 for electrons and the API, AP6, AP7 for high energy protons. All are 
static models which do not consider temporal variations. For the protons this 
is a valid representation because experimental measurements have shown that 
no significant changes with time have occurred. With the exception of the fringe 
areas of the proton belt, that is, at very low altitudes and at the outer edges of 
the trapping region, the possible error introduced by the static approximation 
lies well within the uncertainty factor of 2, attached to the models. Consequently, 
the proton models may be applied to any epoch and no updating process is needed. 

The same is not true for the electrons, which are strongly time dependent, both 
in the inner as well as in the outer zone; especially in the case of the AE2 model, 
which describes the environment as it existed back in 1964, when the artificial 
component of the total population was predominant in the inner zone and solar 
minimum conditions prevailed in the outer zone. Specifically, in 1964 a large 
proportion of the Starfish electrons were still trapped in the belt between 
L = 1.2 and L = 2.0; these seem to have decayed down to natural background- 
levels by 1967-1968. On the other hand, between L = 2.0 and L = 5.5 the 
trapped particle population appears to be strongly solar cycle dependent, in a 
sense that an increase in solar activity towards solar maximum is expected to 
produce an increase in the ambient intensities due to an enhanced trapping 
mechanism. Thus, it is anticipated that the electron fluxes in the outer zone 
may have risen substantially during the years 1969-1970. A practical quantita- 
tive or qualitative treatment of this process is not available yet. In the mean- 
time, estimates of the predicted population-increases vary within the scientific 
community; for a peak somewhere around L - 3.4 ± .2, they range from a 
factor of 30 up to several hundred. 

To obtain reasonable electron predictions for any epoch after 1964, when using 
the AE2, it is necessary to update the model data. In order to account for the 
change in the inner zone, the 1964 fluxes were exponentially decayed to July 1967 
with experimentally measured decay constants. However, lacking suitable data, 
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the electron calculation could not be corrected for solar cycle variations. To 
partially compensate for this, the uncertainty factor attached to the results has 
been raised appropriately. 

Units: 

The results, relating to omnidirectional, vehicle encountered, integral, trapped 
particle fluxes, are presented in graphical and tabular form with the following 
unit convention: 

1. Daily averages: total trajectory integrated flux averaged into: 

particles/cm 2 day 

2. Totals per orbit: non-averaged, single-orbit integrated flux in: 

particles/cm 2 orbit 

3. Peaks per orbit: highest orbit-encountered instantaneous flux in: 

particles/cm 2 sec 

where 1 orbit = 1 revolution. 


Uncertainties : 

We wish to emphasize the fact that the data presented in this report are only 
approximations. We do not believe the results to be any better than a factor of 
2-3 for the protons and a factor of 4-5 for the electrons. It is advisable to 
inform all potential users about this uncertainty in the data. 


Orbit- integrated Fluxes and Spectra: 

The results of our calculations are summarized in Tables 1 (electrons) and 
2 (protons) for all requested orbital configurations. The spectral distribution 
for selected energies is additionally expressed in a percentage form, normalized 
to the E > .5 Mev fluxes for electrons and the E > 5. Mev fluxes for protons. 

A graphical superposition of the data is given, for comparison, in Figures 1 
and 4, and the selected set of integral energies are plotted versus altitude and 
inclination in Figures 2, 3 for electrons and in Figures 5, 6 for protons. 
Classification of orbit-integrated spectra as hard or soft is relative, based on 
an overall evaluation of near earth space in terms of circular trajectories 
between equatorial and polar orbits. 
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On some preliminary graphs discontinuities appeared in the proton spectra. 
These "breaks" occurred because the complete proton environment is being 
described by three (formerly four) independent maps or grids, each valid only 
over a limited energy range; for certain critical orbital configurations the dis- 
continuities are then produced when moving from one energy range to another. 
They are caused, in part, by the exponential energy parameter of the model 
which in many instances had to be extrapolated to make up for lacking data and, 
in part, to insufficient experimental measurements over some areas of B/L- 
space; furthermore, the discontinuities reflect the fact that the available data 
cannot be completely matched at their overlap. In order to overcome such 
spectral breaks, a continuous weighted mean curve was drawn, connecting the 
adjacent segments; it should be regarded as an approximate spectral distri- 
bution. In doing this, the API results [30 < E (Mev) < 50] had to be totally 
ignored sometimes. 

In Figure 1, a rise in the vehicle-encountered fluxes with height may be ob- 
served; this is typical for these altitudes. Also, for a given fixed altitude level, 
a clustering of the inclination curves is evident. A comparison with Figures 
2 and 3 indicates a slight hardening of the spectra for higher altitudes, while 
an increase in inclination apparently produces a small softening effect. Both 
variations are marginal and may be neglected. The electron spectra may be 
classified as moderately hard for near earth space missions. 

In Figure 4, the rise in the accumulated fluxes with altitude shows up again and 
the grouping effect is present. The curves indicate a substantial hardening of 
the spectra with an increase in height. Over the specified range of inclinations 
no spectral change is perceptible, but test calculations with i = 28° for 556 km 
and 370 km produced a sizr able drop of the fluxes in the lower end of the 
spectrum; in addition, the 370 km orbit yielded high-energy flux-values that 
fell below the 52° level. Figures 5 and 6 corroborate these findings. The proton 
spectra may be classified as hard to very hard, in the cases where i - 28°. 

Attached are 24 computer plots depicting the characteristic electron and proton 
spectra of the twelve trajectories, individually. 


Percent of Total Lifetime r spent in the Trapped Particle Radiation Belt: 

The Van Allen radiation belt extends, according to prevailing opinion and 
experimental evidence, to approximately synchronous altitudes along the mag- 
netic equator, that is, to about 6.6 earth radii; there ordered particle trapping 
breaks down. For our purposes, we assume the stable trapping region to be 
bounded and well defined in conventional magnetic B/L-space. As explained 
before, we divide the L-range of the space into an "inner zone" (1.1 < L < 2.0) 
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and into an "outer zone" (L > 2.0). Correspondingly, we will denote the per- 
centage of total lifetime spent by the vehicle in each zone by r* and r°. The 
percent duration spent outside the radiation belt (L > 6.6) will be given by r e 
(r-extemal). For any fixed altitude and inclination r* + t 0 +t' = 100%. 

The calculated values for t‘, t°, and r* are presented in Table 4; they are 
plotted in Figure 9 versus altitude and in Figure 10 versus inclination. 

It is evident from Figure 9 that the r's are not very sensitive to changes in 
altitude; this holds for all tested inclinations. Because of the relatively small 
variations in the r's we may eliminate altitude as a variable and use mean 
values of r's instead, averaged over the altitude range specified for this study. 
These values are plotted in Figure 10. 

As would be expected, the curves show that for low inclinations, i < 30°, the 
vehicle spends its entire lifetime in the inner zone. When inclination is raised 
from 30° to 50° it enters the outer zone for rapidly growing periods of time; 
finally, when inclination exceeds 50°, it begins leaving the trapped radiation 
belt, initially for short intervals of time, which however increase in almost an 
exponential fashion, as i = 90° is approached. In this last stage the satellite 
performs a complete sweep through magnetic L- space, which constitutes the 
transverse motion mentioned earlier, executed twice during each revolution 
(orbit). I 

If additional data were available for the higher inclinations, we believe the 
plotted curves would continue to extend monotonically towards i = 90° with the 
exception of the outer zone contour which indicates a folding at i = 65°. We 
would estimate the final values for a polar orbit to be approximately: 


r i = 45%, t° = 20%, T e = 35%. 


Because the information contained in this section may also prove useful in the 
selection of an orbit-inclination, we wish to bring up and submit for prior con- 
sideration the following related points: 

(a) lasting solar cycle effects, that is, significant increases in the trapped 
electron population from solar minimum to solar maximum, are more 
severely experienced in the outer zone; 

(b) energetic artificial electrons from high altitude nuclear explosions 
(Starfish) have displayed a remarkable longevity, but only in the inner 
zone; there, they contaminated the environment for over 5 years, while 
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they rapidly decayed to background levels in the outer zone (within 
weeks to months). A planned or accidental explosion of another atomic 
device with the appropriate yield and at the right latitude and altitude 
may, very likely, produce conditions similar to those experienced with 
"Starfish", transforming the inner zone again into a radiation hotbed. 


Percent of Total Lifetime r Spent in High Intensity Regions of the Van Allen 
Belt and Percent of Total Daily Flux Accumulated During This Time : 

We arbitrarily define as "High Intensity" those regions of space, where the 
instantaneous, integral, omnidirectional, trapped particle flux is greater than 
10 3 protons with energies E > 5 Mev and greater than 10 s electrons with 
energies E > .5 Mev. 

Table 3 gives the total lifetime percentage r hi for the high intensity regions and 
the respectively accumulated fluxes J, for protons and electrons. These values 
are graphically presented in Figures 7 (t and J vs. i) and 8 (r and J vs. h). 

An interesting feature in Figure 7 is the clustering of the proton curves; 
evidently, the encountered high intensity proton fluxes are almost independent 
of inclination, while the corresponding time-curves show a substantial variation 
in i, which implies that the vehicle passes through regions of successively higher 
absolute intensities as inclination increases, but for shorter total intervals of 
time. This relationship seems to be more pronounced for the upper altitude 
levels. 

A similar features appears also for the electrons at 556 kilometers. This 
independence is clearly expressed in Figure 8 by the zero slope of the altitude 
curves. 


Percent of Total Lifetime r Spent in Flux-Free Regions of Space: 

In the context of this study, the term "Flux-Free" applies to all regions of space 
where trapped- particle intensities are less than one electron or proton per 
square centimeter per second, having energies E > .5 Mev and E > 5 Mev, 
respectively; this includes regions outside the radiation belt. We denote the 
percentage of mission lifetime spent in flux-free regions by r f f . Values of r f f 
for both types of particles are given in Table 5; they are plotted versus altitude 
and inclination in Figures 11 and 12, respectively. 
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It should be noted that the curves for protons and electrons in Figure 11 dis- 
play an inverse relationship of r ff to inclination, that is, as inclination is 
increased the value of r ff decreases for electrons but increases for protons. 
This is consistent with the geomagnetic geometry of the radiation belt; it re- 
flects the fact that the high energy (E > 5 Mev) proton population occupies a 
much smaller volume than the low energy (E > .5 Mev) electron population, 
which in addition displays the so called "horns" (Stassinopoulos, NASA SP-3054) 
at the altitudes under consideration. Consequently, an increase in inclination 
moves the trajectory farther away from the high intensity regions of the proton 
belt into flux-free space, but apparently it still immerses the orbit deeper into 
the electron belt. 

A better picture of this relationship is conveyed by Figure 12. Evidently, as the 
inclination approaches 75° - 80°, the vehicle experiences a maximum exposure 
to accountable electron intensities (r ff getting smaller), whereas the proton 
exposure decreases continuously (r ff becomes larger), gradually levelling off 
to some final limiting value for a polar orbit (probably around r ff % 95%). 

It is reasonable to assume that the electrons will reverse their trend somewhere 
around 76° - 78° and yield possibly larger r ff 's while approaching a polar 
inclination. 

An extension of the curves towards lower inclinations (30° < i < 50°) would most 
likely follow the pattern established in the plot. However, there is reason to 
believe that in the neighborhood of i = 30° a "folding" will occur for the electrons 
as the T f f 's slowly settle to a final limiting value for an equatorial trajectory. 

Obviously, if an inclination were to be selected on the basis of r ff alone, that 
is, total time available in flux-free regions of space, a prior evaluation and 
comparison of the radiation hazards due to the predicted electron and proton 
fluxes would be essential, either in regards to the entire mission or in regards 
to specific important mission requirements. While the proton intensities are on 
the average two orders of magnitude smaller than the electrons, and while they 
afford more flux-free time Tor hi&nci inclinations, their greater mass and harder 
spectra may prove more d maging to t ie mission than the more numerous 
electrons with their lesser flux -free t.me. 

It appears that the dependence of r n on altitude is almost linear for all incli- 
nations and particles. 

Peak Fluxes per Orbit : 

Attached are 24 additional computer plots for the twelve initial trajectories, 
showing the vehicle encountered instantaneous peak electron (E > .5 Mev) and 
proton (E > 5. Mev) intensities per orbit for a sequence of about 30 revolutions. 
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On all graphs a periodic pattern emerges that indicates a daily cycle of about 
15 orbits which may shift slightly, alternating between 14 and 16 orbits; this is 
due to the relative orbit- period, which determines the precession of the 
trajectory. If there were an integer number of orbit-periods in a day, then the 
cycle would be an exact 15 revolutions. For circular flight paths, the orbit- 
period is a function of the altitude. Inclination has no bearing on this process. 

It affects, however, the incident peak fluxes of the electrons and the protons, 
but apparently in different ways, as will be discussed in the next two paragraphs. 

Curves of both types of particles display certain "troughs" or depressions, 
sometimes to the zero- flux level; usually these minima occur at the ends of 
the cycles and may extend over several orbits; occasionally they appear also 
within a cycle, as for instance in the high inclination proton plots. 

In the case of our study, such intra-cyclic flux-less orbits are observed only 
in the proton medium rnd they appear to be a function of the inclination. At 
i = 52° there is but a small depression in the top segment of the contour; at 
i = 65° the trough has already formed, and at i = 72° it is complete. It lasts 
for 1-2 orbits, the most, and seems to hold for all altitudes. 

An increase in height raises the entire curve, thereby lifting the trough. 
Eventually, at some altitude the minimum wili have contracted to within an 
order of magnitude from the top of the curve. 

The inter-cyclic flux-less orbits are apparently a fundamental characteristic 
of the E > 5 Mev proton environment, that is, they show up consistently for all 
altitudes and inclinations. Their number varies from to v 5. The ’’ectrons 
demonstrate such a minimum only at the lowest inclinations considereu, .nd 
only for about 1 orbit. 

Time- and Flux- Histograms: 

Finally, for each of the twelve flight paths, two more plots are includes, one 
for electrons and one for protons, depicting the characteristic averaged 
instantaneous intensities of the trajectory in terms of constant L- bands of 
.1 earth radius width; the percent of total lifetime spent in each L-interval is 
shown on the same graph by the contour marked with the x’s. 
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